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Abstract. In Roscoe 1999a, it was described how the modelling of a small sample of optical rotation curves (ORCs) given 
by Rubin et al 1980 with the power-law Vrot = AR", where where the parameters {A, a) vary between galaxies, raised 



the hypothesis that the parameter A (considered in the form In A) had a preference for certain discrete values. This s pecific 
hypothesis was tested in that pa per ag ainst a sample of 900 spiral galaxy rotation curves measured by Mathewson et al 



1992 



but folded by Persic & Salucci 1995, and was confirmed on this large sample with a conservatively estimated upper bound 



probability of 10^^ against it being a chance effect. In this paper, we begin by reviewing the earlier work, and then describe 
the an alvses of three additional samples; the first of these, of 1200+ Southern sky ORCs, was published by Mathewson & Ford 
1996, the second, of 497 Northern sky ORCs, is a compos i te sam ple provided by k ind pe rmission of Giovanelli & Haynes 



published in the sequence of papers D ale et al 1997, 1998, 1999 and Dale & Uson 2000|, whilst the third, of 305 Northern 



sky ORCs, was published by Courteau 1997. These analyses provide overwhelmingly compelling confirmation of what was 
already a powerful result. Apart from other considerations, the results lead directly to what can be described as a 'second 
generation of Tully-Fisher methods.' We give a brief discussion of the further implications of the result. 
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1 . Introduction 

This paper describes the analyses of four large optical rotation 
curve samples to show how the hypothesis that 'spiral galaxies 
are constrained to occupy discrete dynamical classes' is sup- 
ported by the data as a statistical certainty. At a practical level, 
the result has considerable immediate implications for the ab- 
solute determination of zero points for classical Tully-Fisher 
methods and, in the wider context of the overall analysis, pro- 
vides a class of second generation Tully-Fisher methods in 
which absolute Tully-Fisher calibrations and linewidth deter- 
minations can, in principle, be simultaneously determined on 
any given sample of ORCs. At a deeper level, the basic result 
appears to have implications for our understanding of galaxy 
evolution. The result is so unexpected, that a short review of 



already published material (Roscoe 1999a) is likely to be use- 
ful to the present reader 

Arguments based on certain symmetry considerations - the 
nature of which are not immediately relevant here - led us 
to consider the possibility that the segments of optical rota- 
tion curves which occupy disc regions (given an operational 
definition later in this text) might be reasonably described - 
in an overall statistical sense - by power laws in the form 
Vrot — AR"', with the parameters A and a being deter- 
mined empirically for each galaxy in turn. As a means of 
gaining familiarity with this idea, we considered the small 
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sample of 21 ORCs published by Rubin et al 198C from this 
point of view. Of this sample of 21 ORCs, only twelve man- 
ifested reasonably monotonic behaviour and so were selected 
on these grounds alone as reasonable candidates for a power 
law analysis. Subsequently, a linear regression of the model 
In Vrot = In ^ + a In i? onto each of the twelve ORCs pro- 
vided twelve sets of parameter-pairs (a, In^). The first clear 
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Table 1. 


Twelve Rubin, 


Ford & Thonnard 1980 








Galaxy 


In A 


Galaxy 


\nA 


N3672 


3.6 


U3691 


3.6 


N3495 


4.0 


N4605 


4.0 


10467 


4.1 


N0701 


4.1 


N1035 


4.1 


N4062 


4.5 


N2742 


4.5 


N4682 


4.5 


N7541 


4.6 


N4321 


4.9 




Table 2. 


In A data 




RFT 


Pred value Actual value 


scale 


with MFB 


MFB 




scale 




scale 




3.5 


3.81 




3.85 


4.0 


4.22 




4.24 


4.5 


4.63 




4.72 


5.0 


5.04 




5.06 





result of this mini-analysis was that a and In A appeared to be 
very strongly correlated - and this particular as pect h as now 
been analysed in detail using Persic & Salucci's |l995| folding 
solution for 900 ORCs from the Mathewson et al 1992 sample 
(Roscoe [l999b| ). 

However, as reference to table |l| shows (the entries of 
which have been rounded to the nearest decimal), a curi- 
ous numerical coincidence arose - specifically, that every one 
of the twelve In A values lay between ±0.15 of an integer 
or half-integer value - a coincidence that has odds around 
1 :500 of being a chance occurrence ( a-posteriori probabil- 
ities !). Of course, the integer/half-integer values themselves 



can be of no possible significance since, if Rubin et al 1980 



had estimated distance scales using a value of H significantly 
different from the 50km/sec/Mpc they actually used, then a 
completely different set of In A values would have resulted. 
So, the coincidence was simply that of regularity in spac- 
ing which would probably have not been noticed with, say, 
H = 70km/sec/Mpc. Anyway, curiosity provided a suffi- 
cient motivation to consider the matter further, using the Persic 
& Salucci |1995| sample of 900 folded ORCs. This sample 
had its distance-scaling determined by a Tully-Fisher relation- 
ship calibrated by Mathewson et al 1992 against Fornax using 
H = 85km/sec/Mpc, so that the integer/half-integer hypoth- 
esis for In A is not appropriate. However, a simple analysis 



(described in Appendix B of Roscoe 1999a, and relying on 



the inv estigation of the {a, In A) coiTelation given in Roscoe 



1999t ) reveals the relation 



IuAmfb ~ 0.82 IuArft + 0.94 

where Amfb denotes the value of A determined using the 
Mathewson et al 1992 scaling, whilst Aupx denotes its value 
determined using the Rubin et al 1980 scaling. Using this lat- 
ter relation, the integer/half-integer values of In A in the Rubin 
et al 1980 scaling transform into their corresponding value in 



the Mathewson et al 1992 scaling according to the first two 
columns of table 0. The actual In A distribution of the 900 
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Fig. 1. In A distribution for the Mathewson et al 1992 sample 



with Persic & Salucci folding. Vertical dotted lines indicate 
actual peak centres. Vertical solid bars indicate predicted peak 
centres from the Rubin et al 198C sample. Bin width = 0.055. 



ORCs of the Mathewson et al 1992 sample folded by Persic & 
Salucci 1995 is given in figure 111. The vertical solid bars indi- 
cate the predicted positions of the peaks, based on our analysis 
of the Rubin data, as in the second column of table 0, whilst 
the vertical dotted lines indicate actual peak centres, as in the 
third column of table 0. The correspondence between the peak 



positions, predicted on the basis of the twelve Rubin et al 198C 
galaxies of table 111 and the actual peak positions is clearly re- 
markable. A crude, but extremely conservative, upper bound 
estimate of the probability of the peaks in the distribution of 
figure H occurring by chance, gi ven th e original hypothesis 
defined on the small Rubin et al 1980 sample, was given in 
Roscoe |l999a| as 10"'^. 



The implications of this result are so potentially signifi- 
cant, that it has become essential to test the specific hypothesis 
against new samples. This is done using three additional sam- 
ples in the following sections, and the results are overwhelm- 
ingly in favour of the hypothesis. That is, it very much ap- 
pears as though we are seeing evidence for discrete dynamical 
classes in spiral galaxy discs. 

1.1. Organization of paper 

The details of the samples used are described in §H whilst es- 
sential computational details are described in §|| Sections §H 
§0, §g| and §^ describe the core analyses of the four samples, 
whilst the statistical analysis of the results of these analyses 
is described in §|] We discuss the implications of the analysis 
for Tully-Fisher methods in §^ Major stability issues are dis- 
cussed in §[lQ whilst a brief discussion of possible theoretical 
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implications is given in mU The whole is summarized in §0 
Various secondary issues are covered in the appendices, and 
referred to as necessary. 

2. The Samples 

The basic relevant characteristics of the four samples anal- 
ysed are given in Table g, listed in order of probable quality 
as judged by either mean apparent magnitude, or by % of late- 
type spirals (which have a higher hydrogen content than early- 
type spirals, and are therefore likely to be associated with more 
accurate Ha Doppler shift measurements). We discuss, and 
analyse, the samples in order of likely quality. 

There is an additional potential problem pointed out 
by Bosma 1992 in the particular context of the sample of 



DGHU. The rotation curve measurements were done on the Mt 
Palomar 5m telescope and the CTIO 4m telescope, whilst the 
/-band photometry was done on the KPNO and CTIO 0.9m. 
telescopes. 

Although this sample is, typically, three to four times more 
distant that the objects of the first sample the ORC determina- 
tions were made with much larger telescopes, and the sample 
has a very high proportion of late-type objects. We can there- 
fore expect that the quaUty of the sample will bear comparison 
with both of the Mathewson et al samples. We note that about 
11 % of this sample consists of edge on galaxies. 



2.3. The Third Sample, Mathewson & Ford \1 996 



Mathewson et al 1992. Specifically, he points out that about 



17% of the Mathewson et al galaxies are seen edge-on and 
that, for these galaxies, internal absorption effects mean that 
ORC measurements are subject to relatively large errors. Table 
lists the % of edge-on galaxies in each of the samples. We 
discuss the effects of excluding these edge on galaxies when 
the detailed analyses of the samples are given. 



The thi rd sam ple of 1200h- ORCs was obtained by Mathewson 
& Ford 1996 in the period 1991-93 as part of the same observ- 
ing programme that gave the original 965 ORCs of Mathewson 
et al 1992. The main differences between the Mathewson & 



2. 1. The Original Sample, Mathewson et al 1992 



In the period 1988-90, Mathewson et al |1992| measured Ha 
and Nji rotation curves for 965 Southern sky spirals on the 
2.3m telescope at Siding Spring Observatory, whilst the cor- 
responding /-band photometry was obtained using the Im and 
3.9m telescopes. The Njj observations were used to provide 
an estimate of the internal measurement accuracy of the //„ 
observations, and these estimates, in the form of a parameter 
varying on the range (0, 1), were provided for each velocity 
measurement on every ORC. So far as we are aware, the Njj 
observations have not been made available. 

Persic & Salucci took this sample of 965 ORCs and sub- 
sequently produced a sample of 900 good-to-excellent quality 



For d|1996| and Mathewson et al |1992| samples are given in ta 
ble |: It is clear that the Mathewson & Ford 1996 sample is, 
on average, 73% more distant than the Mathewson et al 1992 
sample, meaning that, on average, we only receive 1/3 as much 
light (all other things being equal) from each of the objects; 
this is consistent with the fact that there is an average of a 
1.1 apparent magnitude difference between the samples. This 
large difference in 'light received' indicates that we can expect 



ORC measurements on the Mathewson & Ford 1996 sample 



to be significantly less accurate than those on the Mathewson 



etal 1992. 



Furthermore, the Mathewson & Ford 1996 sample consists 
of only 18% late-type spirals compared to 43%, 57% and 45% 
respectively for the other three samples. Since late-type spirals 
are significantly richer in hydrogen than are early-type spirals, 
and since ORCs are measured primarily in Ha, then we can 
expect the quality of velocity measurements in the Mathewson 



& Ford 1996 sample to rank behind that of the other samples 



folded ORCs (Persic & Salucci 1995), suitable for their pur 



for this reason also. We note that only about 6% of this sample 
consists of edge on galaxies. 



pose of modelling the internal dynamics of spiral galaxies. It 
was on this sample that the 'discrete dynamical classes' hy- 



pothesis was originally tested (Roscoe 1999a), and which is 
represented in figure 111 We note that about 17% of the sample 
consists of edge-on galaxies. 

2.2. The Second Sample, Dale, Uson, Giovanelli & 
Haynes \1997\ et seq 



2.4. The Fourth Sample, Courteau 1997 



The second sample of 497 ORCs was provided by the kind 
permission Riccardo Giovanelli and Martha Haynes specifi- 
cally for the present analysis, and is a composite of samples 



The fourth sample, of 305 ORCs, was selected by Courteau 
from a sample of Sb, Sc field galaxy ORCs (Courteau 1997) 
for a linewidth/Tully-Fisher study, and differs from the first 
three samples in being the only sample with /i'-band photome- 
try. The original observations were made using the Shane 3m. 
telescope at Lick Observatories and the du Pont 2.5m tele- 
scope at Las Palmas. 

As reference to table ^ shows, the Courteau 1997 sample 
is almost as distant as is the Mathewson & Ford 1996 sam- 



pubUshed by Dale, Giovanelli & Haynes 1997, 1998, 1999 



and Dale & Uson 2000, and was originally selected for stud- 
ies of peculiar velocity distributions. As occasion demands we 
shall refer to the sample as either that of Dale et al, or that of 



pie, but it contains a similar proportion of late-type spirals to 
that contained in Mathewson et al 1992. Also, the telescopes 
used were a little larger and so, all other things being equal, 
we would expect the quality of these ORCs to be midway be- 
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Table 3. Comparison of the four samples 







Mean 


Mean 


% Late 


% Edge 




Sample 


Sample 


distance 


apparent 


type 


on 


Telescope 




size 


km/sec 


magnitude 


spirals 


discs 


diameter 



MFB 1992 


900 


3651 


12.3 (J-band) 


43% 


17% 


2.3m 


DGHU 1997 


497 


13747 


15.0 (J-band) 


57% 


11% 


4m, 5m 


MF 1996 


1200 


6311 


13.4 (/-band) 


18% 


6% 


2.3m 


SC 1997 


305 


5854 


13.5 (i?-band) 


45% 


<1% 


2.5m, 3m 



tween that of Mathewson et al 1992 and Mathewson & Ford 



1996 and similar overall to the Dale et al sample. However, the 



sample is considerably smaller than the others so that statistics 
performed upon it will be correspondingly less significant. We 
note that less than 1 % of this sample consists of edge on galax- 
ies. 

3. Miscelleneous computational details 

3. 1. Pre-folding data filter 



Persic & Salucci 1995 found that an ORC could only be folded 
with sufficient accuracy for their purpose if individual veloc- 
ity measurements not satisfying a pre-determined accuracy 
condition were rejected. This is also our experience with the 



auto-folder method described in Roscoe 1999c. Accordingly, 



an individual velocity measurement on any given ORC is re- 
tained only if it has an estimated absolute error < 5%. For the 



Mathewson et al 1992, Mathewson & Ford 1996 Dale et al 



1997 et seq and Courteau 1997 samples this process leads to 
losses of 35%, 25%, 46% and 46% respectively of all individ- 
ual velocity measurements. The folding process for each ORC 
is performed once this data-filtering process is completed. But 
the data-filtering itself inevitably means that some ORCs are 
left with insufficient velocity points to permit subsequent ac- 
curate folding. The attrition rates for ORCs lost to the overall 
analysis via this process are given by 3%, 4%, 8% and 7% for 
each of the four samples respectively. 

3.2. The computation of In A 

The 'discrete dynamical classes' hypothesis is a statement 
which specifically concerns the values assumed by the set of 
\nA parameters, computed for each folded ORC in turn. It 
is therefore necessary to state clearly how this parameter is 
computed. The following discussion assumes that each ORC 
is folded and corrected for inclination as a matter of course. 

In the original mini-analysis of the 21 ORCs of Rubin 
et al (1980), we rejected nine on the grounds that they were 
strongly non-monotonic, and therefore not amenable to a 
power-law analysis. However, the application of any such sub- 
jective culling procedure to the present large-scale analysis 
would seriously compromise its objectivity, and so some algo- 
rithmically defined 'black-box' technique is required to per- 
form the equivalent task. 



Rather than trying to identify complete ORCs which are 
unsuitable in some way for our analysis, we note that most 
non-monotonic behaviour is on the interior parts of ORCs 
and, accordingly, develop an algorithmic statistical technique 
to identify 'bad' interior sections, say < R < i?mm, where 
they exist. The whole of these 'bad' regions are then rejected 
and the whole remaining ORC is used in our analysis without 
further processing. Apart from the ORC attrition arising from 
the pre-folding data-filter, the four ORC samples are used in 
their entirety. 

3.3. Thie algorithim 

The process to be described uses the techniques of linear re- 
gression and, following conventional definitions, an observa- 
tion is reckoned to be unusual if the predictor is unusual, or if 
the response is unusual. For a p-parameter model, a predictor 
is commonly defined to be unusual if its leverage > 3p/N, 
when there are N observations. In the present case, we have 
a two-parameter model so that p = 2. Similarly, the response 
is commonly defined to be unusual if its standardized residual 
> 2. The computation of {a, In A), for any given folded and 
inclination-corrected rotation curve, can now be described as 
follows: 

1. Form an estimate of the parameter-pair (a, In A) by re- 
gressing the In Vrot data on the In R data for the folded 
ORC; 

2. Determine if the innermost observation only is an unusual 
observation in the sense defined above; 

3. If the innermost observation is unusual, then exclude it 
from the computation and repeat the process from (1) 
above on the reduced data-set; 

4. If the innermost observation is not unusual, then no further 
computation is required - the current values of (a, In A) 
are considered as final. 

This algorithm has the result that, on average, {a,\nA) is 
computed on the exteri or 88% of the data points in each ORC 
of the Mathewson et al 1992 sample, the exterio r 87% of data 
points in each ORC of the Mathewson & Ford 1996 sample, 
the exterior 91% of data points in each ORC of the Dale et al 
1997 et seq sample and the exterior 91% of the data points in 



each ORC of the Courteau 1997 sample 
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Whilst this process has an obvious statistical objectivity, 
it is not clear that it has any basis in physics. However, using 
Rmin to denote the radial position of the innermost remain- 
ing point of any given ORC after the application of this pro- 
cess, it turns out that Rmin is a very strong, but noisy, tracer 
of the optical radius, Ropt, as defined, for example, by Persic 
& Salucci 1995 The implication of these correlations, given 
their strength, is that the computed Rmin has some physi- 
cal significance - possibly as a transition boundary between 
bulge-dominated and disc-dominated dynamics. The details of 
these correlations, for all the samples analysed, are given in 
Appendix y. 

3.4. The representation of In A 

All the In A frequency diagrams shown in this paper are ob- 
tained using the same bin-width (A In yl — 0.055) and ini- 
tial point (In A = 2.2) that were used in the original paper, 
Roscoe 1999a, on this topic. There are therefore no hidden 
degrees of freedom available to enhance the signals being dis- 
cussed. 



3.5. The absolute determination of TF zero points 

Tully-Fisher methods are central to our analysis and, as is well 
understood, whilst TF calibrations provide an absolute deter- 
mination of the gradient term, the value of the zero point is 
considered to depend on the value of H assumed for the gen- 
eral calibration procedure. 

In accordance with this, each sample is analysed by con- 
straining our TF gradients to lie within the quoted error bars 
of the absolute TF gradient determinations made by the as- 
tronomers who compiled each of the various samples. We then 
consider if, for each sample, a zero point can be determined 
which gives rise to the significant peak structure of the hypoth- 
esis. If the answer to this question is positive then, in practice, 
a method for the absolute determination of the zero points will 
have been demonstrated. 



3.6. A simple diagnostic guide for TF calibrations 

For ideal data, for which no systematic bias of any kind ex- 
ists, we should find a statistical equality between Hubble mag- 
nitudes and Tully-Fisher magnitudes; that is, we should find 
AfrpF ~ Mnubbie over the magnitude range of the sample 
whenever a sample is quiet in the Hubble sense. On the basis 
of this latter assumption, we apply the following simple diag- 
nostic guide as a means of making quick - and usually reliable 
-judgements in our various analyses: 

1 . Compute the regression model Mtf = AM Hubble + -B on 
the sample for some assumed value of H\ 

2. Adopt some standard reference range on Hubble magni- 
tudes. For example, we use (—23.3,-18.2) which con- 



tains about 95% of the Mathewson et al 1992 Hubble mag- 
nitudes with H = 85km/sec/Mpc. 
3. Use the regression model to compute the magnitude map- 
ping 

( — 23.3, — 18.2)/f„f,;,ie — > {Mmin,Mmax)TF, 

after rejecting outliers, and use this mapping to make qual- 
ified judgements about the TF calibration relative to the 
assumed value of H. 



3. 7. Other Essential Routine Issues 

There are two other essential, but routine, issues that are dealt 
with in the appendices. These are: 

- Tully-Fisher methods are fundamental to this analysis and 
so it is necessary to certain that Tully-Fisher scatter cannot 
wash out the signals being claimed. This is dealt with in 
appendix §^. 

- There is always the possibility that the hypothesised effect 
is an artifact. The various possibilities include the origi- 
nal process of measuring ORCs, the methods of linewidth 
estimation and methods of folding. These possibiUties are 
discussed in §g} 



4. The Analysis of the Mathewson et al |1992 
Sample 



The M athewson et al |1992| sample, like the Mathewson & Ford 



1996 sample, is drawn from an area of the sky which Lynden- 



Bell et al |1989| believe contain the Great Atti-actor (GA) and 
approximately one half of the Mathewson et al 1992 and 
Mathewson & Ford 1996 samples lie within (Mathewson et al 
1992| definition of) the GA region, 260° < I < 360° -40° < 
< 45°. In their figure 12, Mathewson et al 1992 use their 



Tully-Fisher calibration to show that inside the GA region, the 
data exhibits a clear bias consistent with some form of large- 
scale flow whilst, outside of the GA region, there is no such 
bias. For this reason, we restrict applications of the magnitude 
mapping diagnostic on this sample to the non-GA regions. 

4. 1. The Mathewson et al calibration for MFB data 



Mathewson et al 1992 calibrated their Tully-Fisher relation 
against the Fornax cluster (for which there is a very narrow 
redshift dispersion) on the basis of the assumption that Fornax 
is at 1340km/sec (using H = 85km/sec/Mpc). Uniquely 
amongst the samples analysed here, their linewidth determi- 
nations were made using an intuitive case-by-case 'eye-balF 
method (private communication). They obtained 



M = -8.18 (±0.12) log V;, 



(1) 



as their Malmquist bias corrected TF form. In fact, Mathewson 
et al do not report error bars on their TF gradient, but do re- 
port error bars of ±0.26 mags for magnitude determinations in 
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Fig. 2. In j4 distribution for the Mathewson et al 1992 sample 
with auto-folding and original Mathewson et al TF scaling; 
Vertical dotted lines indicate peak centres of Persic & Salucci 
solution. Bin width = 0.055 
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Fig. 3. In j4 distribution for the Mathewson et al 1992 sample 
with auto-folding and rescaled TF; Vertical dotted Unes indi- 
cate peak centres of Persic & Salucci solution. Bin width = 
0.055 



Fornax. Our estimate of ±0.12 as the error bar for their gradi- 
ent determinations is derived from this. Using their mean value 
gradient of —8.18 and the same zero point, with Mathewson's 
assumed H = 85 km/sec/Mpc, we find that, for the subsample 
exterior to the GA region, the diagnostic magnitude mapping 
gives 



(-23.3,-18.2)^85 ^ (-23.0,-18.1)tf 



(2) 



which indicates a possible discrepancy at the bright end be- 
tween Hubble and Tully-Fisher magnitudes in the non-GA re- 
gion. 

The \\iA frequency distribution for the Mathewson et al 
sample with the calibration dl]) (with Mathewson's quoted 
—8.18) generated by our own folding technique is shown in 
figure 0. The vertical dotted lines in this latter figure mark 
the positions of the peaks A, B, C, Z? of the Persic & Salucci 
1995 solution, in figure [l|. It can therefore be seen that the 



peak structure revealed by the Persic & Salucci |1995| folding 
method is not an artifact of th eir me thod, but is an objective 
feature on the Mathewson et al 1992 sample. 



4.2. The effect of rescaling the TF relation within 
the error bars 

According to the diagnostic magnitude mapping, (g), the ref- 
erence range of TF magnitudes generated by the scaling (|l]) is 
compressed and shifted to the dim end relative to the Hubble 
magnitudes. This suggests that a possible rescaling designed 
to expand the reference range of TF magnitudes could be con- 
sidered. Such an expansion, consistent with Mathewson et al's 
original calibration, is obtained by steepening the TF gradient 
within the quoted error bars. We find that, if the peak struc- 
ture is to be optimized, then a steepening of the gradient (in 



the negative sense) must be accompanied by a reduction in the 
zero point. The overall process is completed by an 'eye-ball' 
iteration which can be described as follows: 



Guess a new gradient value within the quoted error bar; 
Adjust the zero point to maximise the peak structure for 
this gradient; 

Perform the diagnostic magnitude mapping, adjusting H 
as necessary to get the best match at any given iteration; 
Repeat until a satisfactory magnitude mapping is obtained. 



4, 
We find 

M = -8.30 log V^rot -2.56, 
gives the magnitude mapping 

(-23.3,-18.2)^95 ^ (-23.2,-18.2)tf 



(3) 



(4) 



computed with H — 95 km/sec/Mpc. There is now no signif- 
icant magnitude bias or shifting. The corresponding In A fre- 
quency map is given in figure g, and we see a small increase 
in signal strength on all four peaks. In summary, the combined 
effects of the improved magnitude mapping, together with the 
increase in the strength of the In A peak structure suggest that 
the recalibration (Bj) is justified. 

4.3. The effect of excluding the edge-on galaxies 

In practice, we find that excluding the edge-on galaxies, which 
are about 17% of the sample, has the effect of producing only 
a pro-rata reduction in the heights of the peaks in figure |[ 
In other words, the inclusion of the edge-on galaxies in the 
present analysis has a neutral effect upon it. The reason for this 
result is almost certainly that the effects of internal absorption 
are considerably greater on the interior parts of ORCs than 
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they are on the exterior parts, and it is precisely these parts 
which are discarded in the data reduction process described in 
9a for the computation of (a, In A). 

5. The analysis of the Dale, Giovanelll, Haynes 
& Uson sample 

From our point of view, the value of the work of Dale et al 



which is completely independent in all of its aspects of the 
Mathewson et al samples. 

As a general comment, it is to be noted that part of the 
explicit general programme of Giovanelli, Haynes et al has 
been to develop a highly accurate TuUy-Fisher /-band 'tem- 
plate' relation for the purpose of establishing a cluster inertial 
frame out to z w 0.06. Thus, Giovanelli et al 1997 are able 



to quote the very tight error bars of ±0.02 for their TF gradi- 
ent determinations. To ensure such accuracy, much effort has 
been put into establishing a reliable linewidth estimation tech- 



nique. This technique, which is described by Dale et al |1998 
is an algorithmic approach based on the estimate Vgo% — ^io% 



first introduced by Dressier & Faber 1990 However, it should 
be noted that this definition is most reliable when ORCs ex- 
tend out to at least Rs3, the optical radius as defined by 



Persic & Salucci 1995 . The reason is that, as Dale et al 1998 
point out, the definition then recovers Vopt in a reliable way 



and, as Persic & Salucci 1991, 1995 have shown, Vopt pro- 
vides a reasonable basis for linewidth definitions. To deal with 
those ORCs which were insufficiently extensive. Dale et al 



1998 introduced an elaborate procedure which involved the 



use of Nil measurements to calculate a parameter, called a 
'shape factor', for each ORC. This shape factor was subse- 
quently used to extrapolate ORCs out to Rg^ and to correct 
raw VJ3o% — V^io% hnewidths as necessary. 

5. 1. The Dale et al calibration of Tully-Fisher 



The Dale et al 1998 calibration of the Tully-Fisher relation is 
given by 



Mtf = -7.68 (±0.02) log Kot -4.12 + 5 log/i 



(5) 



where H = 100ft- km/sec/Mpc is undetermined. The error bars 
on the gradient are so tight that, effectively, there is no freedom 
to vary the estimated TF gradient. However, the analysis of the 
Mathewson et al sample in §^ strongly suggest a value H = 
95 km/sec/Mpc so that h — 0.95 in (||). Thus the suggested 
calibration is given by 



Mtf = -7.68 log Kot - 4.23, 



(6) 



and the corresponding In A frequency diagram is given in fig- 
ure Q. We see that the predicted peak structure is strongly 
reproduced. Similarly, the diagnostic magnitude mapping is 
given by 

(-23.3, -18.2)h95 -^ (-23.2, -18.1)tf 



Dale et al sample of 454 objects 
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Fig. 4. The In A distributions for the Dale et al sample, with 
auto-folder and using the Dale et al TF calibration for H = 
95 km/sec/Mpc. Vertical dotted lines indicate peak centres 
of figure |. Bin width = 0.055 



which, apart from a very slight uniform shift to the dim, indi- 
cates the complete absence of any systematic magnitude bias. 



5.2. The effect of excluding the edge-on galaxies 

About 1 1 % of the Dale et al sample consists of edge-on galax- 
ies, and the effect of excluding these galaxies is the same as for 
the previous two samples - that is, there is merely a pro-rata 
reduction in the peak heights. 



6. The Analysis of the Mathewson & Ford |1996 
sample 



The Mathewson & Ford 
al 



1996 sample, like the Mathewson et 



1992 sample, is also drawn from an area of the sky which 



contains the, so-call ed, GA and approximately one half of the 
Mathewson & Ford 1996 sample lies within the GA region, 
260° < I < 360°, -40° < 6 < 45° but, as reference to table | 
sho ws, is a n average of 70% more distant than the Mathewson 
et al 1992 sample and is therefore considerably less bright. 



6. 1. The Mathewson et al calibration for MF data 



A direct application of the Mathewson et al 1992 TF calibra- 
tion, given at (|]), to the Mathewson & Ford 1996 sample gives 
the In A frequency diagram of figure |} This figure shows that, 
whilst the B, C and D peaks are reproduced, the distribution 
is extremely noisy and must be considered poor Given the 
quality of the results obtained from the Mathewson et al 1992 
sample, and shown in figures g and g, we can suppose that 
some kind of problem exists. 
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MF 1996: Original TF calibration 
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Fig. 5. The In A distributions for the Mathewson & Ford 1996 
sample with auto-folding using Mathewson linewidth esti- 
mates and original TF calibration. Vertical dotted lines indi- 
cate peak centres of figure 0. Bin width = 0.055 

As a possible means of understanding what the prob- 
lem might be, we consider the subsample of Mathewson & 



Ford 1996 data which is exterior to the GA region (dis- 
cussed in § lb and calculate how the reference range of Hubble 
magnitudes maps into Tully-Fisher magnitudes for this data. 
Using the original Mathewson TF calibration (for which H ~ 
85km/sec/Mpc), we find the diagnostic magnitude mapping 
for non-GA objects 



-23.3,-18.2)^85 ^ (-23.0,-17.4)tf 



(7) 



which implies the existence of a very strong systematic bias 
which shifts the whole Mathewson & Ford 1996 sample to 
lower luminosities. This is improved at the bright end when 
our modified Mathewson TF calibration (^ is used, but the 
dim-end bias remains essentially unchanged.. Given that the 
non-GA objects are not believed to be participating in any 
large-scale flow, there are two basic possibihties for explaining 
the mapping (Q) which can be listed as 

- The MF Hubble luminosities are very much overestimated 
at the dim end; 

- The MF Tully-Fisher luminosities are very much underes- 
timated at the dim end. 

The fi rst possibility seems unlikely since Mathewson & Ford 



1996| photometry is in the / band for which the internal and 
external extinction mechanisms are well understood, and for 
which well-tested correction techniques exist and have been 
applied by Mathewson & Ford 1996. The second possibility 
would necessarily have its source in the systematic underesti- 
mation of dim-end optical linewid ths for this relatively distant 
sample. Since Mathewson & Ford 1996 (and Mathewson et al 



1992) used an intuitive 'eyeball' technique for linewidth esti- 



mation, this latter possibility seems the most likely explana- 
tion. On the basis of the working hypothesis that Mathewson 



& Ford's linewdith estimates are subject to a systematic bias 
relative to the Mathewson et al 1992 linewidths, there are now 



two possible ways to proceed: 

1 . We can either derive our own linewidth estimates directly 
from the sample using some algorithmic technique; 

2. Or we can take advantage of the idea that, where system- 
atic linewidth bias exists, it can accounted for by a com- 
pensating recaUbration of the Tully-Fisher relation. 

We consider both of these approaches in the following sec- 
tions. 

6.2. TF based on Vopt linewidths 

The first approach is based on the idea of generating our 
own linewidth estimates. The most simple algorithmic tech- 
nique for linewidth estimation is that based on Vgo% — Viq%. 
However, as discussed in §q|, linewidths based on this defini- 
tion can only reliably be used when ORCs extend out to at least 
to the optical radius, and Dale et al used an elaborate proce- 
dure to negotiate this problem for those ORCs which did not. 
This approach is not possible in the present case because the 
required Nu measurements are not available for Mathewson 
data. But Dale et al 



1998 



also report that their definition of 
linewidth recovers Persic & Salucci's Vopt whenever ORCs 
are sufficiently extensive. This implies that using Vopt as a ten- 
tative linewidth definition for Mathewson & Ford data should 
allow a TF calibration very similar to that of Dale et al's given 
at(|). 

However, although R^^ (angular measure) estimates are 
provided with the Mathewson & Ford 1996 sample, the corre- 
sponding Vopt are not. We circumvent this problem by using 
the power-law V = A i?" fitted to the folded, but as yet un- 
sealed, ORC to estimate Vopt at i?83 (angular measure). This 
process is possible prior to scaling simply because a values 
are independent of scaling. With these linewidths, we find that 
the exact Dale et al calibration (pt) with a zero point of —4.08, 

i\/ = -7.681ogl4ot-4.08, 

gives the In A frequency diagram of figure g[ By contrast with 
figure I], we see that the whole A, B, C, D peak structure is 
now very well reproduced. Interestingly, and as judged by the 
absence of any bias in magnitude mapping 

(-23.3,-18.2)^85 ^ (-23.2,-18.3)tf, 

this latter TF calibration corresponds more closely to H = 



85 km/sec/Mpc for the non-GA Mathewson & Ford 1996 sub- 
sample than the higher value of iJ = 95 km/sec/Mpc. 

6.3. Linear reseating of TF for linewidtfi 
bias-correction 

The relative success of the Vopt linewidth estimate provides 
further circumstantial evidence for the working hypothesis 
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Fig. 6. The In A distributions for the Mathewson & Ford 1996 




Fig. 7. The In A distributions for the Mathewson & Ford 1996 



sample with auto-folding using Vopt linewidths and Dale et 



sample with auto-folding using the TF recalibration with lin- 
al TF calibration with shifted zero point. Vertical dotted lines ear rescaling. Vertical dotted lines indicate peak centres of fig- 
indicate peak centres of figure Bin width = 0.055 ure 0. Bin width = 0.055 



that the Mathewson & Ford linewidths are, in fact, systemat- 
ically biased. We make the simplest possible assumption that 
such a bias can be accounted for by a linear recalibration of the 
TF relation designed to ensure that the diagnostic magnitude 
mapping is bias-free. We find that the rescaled TF relation 



MTF = -7A6\ogVrot-^.85, 



(8) 



gives the H — 88 km/sec/Mpc diagnostic magnitude mapping 

(-23.3,-18.2)^88 -> (-23.3,-18.2)tf, (9) 

for non-GA ORCs, which is perfect. Corresponding to this im- 
provement, we find that the recalibration (pb gives the In A fre- 
quency diagram of figure ^ which is likewise a considerable 
qualitative improvement over figure p] In particular, the A and 
B peaks are now perfectly and strongly reproduced, whilst the 
C peak is clearly present, but noisy. The D peak is virtually 
non-existent. 

6.4. The effect of excluding the edge-on galaxies 

Only about 6% of the Mathewson & Ford (1995) sample con- 
sists of edge-on galaxies, and the effect of excluding these 
galaxies is the same as for the Mathewson et al 1992 sample - 
that is, there is merely a pro-rata reduction in the peak heights. 



7. The analysis of the Courteau [1997| sample 



As with the Dale et al sample, the value of the Courteau sample 
lies in its complete independence in all of its aspects of the 
Mathewson et al samples. 

The Courteau 1997 analysis was primarily designed to ad- 
dress the problem of linewidth definitions, with a view to ob- 
taining a standardized objectively defined algorithmic defini- 
tion. Courteau considered several possibilities for linewidth 



definitions, and we present results using his Vmax and V2.2 
definitions (his estimated 'worst' and 'best' respectively). 
Because his explicit concern was to make a comparative inves- 
tigation of several linewidth estimation techniques, less effort 
was expended in establishing particularly accurate TF caUbra- 
tions for each of the linewidth estimates considered. For exam- 



ple, unlike Mathewson et at 1992 or Dale et al 1998 , Courteau 
did not calibrate his Tully-Fisher relations on specially cho- 
sen low-redshift dispersion clusters, but on his whole sample, 
which has a redshift dispersion in excess of 10000 km/sec. It 
follows that his calibrations are n ot likely to be as accurate 
as those of Mathewson et al 1992 and Giovanelli et al 1997 



Thus, Courteau's quoted error bars on gradient determinations 
for each of the linewidth estimators considered are, typi cally, 
±0.2, which is to be compared with Giovanelli et al's 1997 
quoted error bar of ±0.02. 

7. 1. The Courteau calibrations of Tully-Fisher for 
two linewidths 

7.1 .1 . The Vmnr linewidth TF calibration 



The Courteau 1997 calibration of the Tully-Fisher relation for 
his Vmax linewidths, using his H = 70 km/sec/Mpc, is given 

by 

Mtf - -6.09 (±0.26) logVmax - 7.22 

where the gradient error bar is as given by Courteau. We sub- 
sequently find that the recalibration 

Mtf --6. 19 logy™, -7.50 

gives the In A frequency diagram of figure ||. Except for the A 
peak which is non-existent because of lack of data at the dim 
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end of the Courteau sample, we see that the peak structure of 
the hypothesis is very well reproduced - although it is much 
noisier than the peak structures of Mathewson et al and Dale 
et al shown in figures H and H respectively. Note that the recal- 
ibrated gradient is comfortably within Courteau's quoted error 
bar. 

It is also of interest to note that the net effect of the recal- 
ibration is to increase the luminosity of the average galaxy by 
about 0.48 magnitudes. This is well within Courteau's quoted 
error bar for Vmax calibrations of 0.55 magnitudes and so 
there is no evidence to suggest that the recalibration is sig- 
nificantly different from Courteau's original Vmax calibration 
which assumed H = 70km/sec/Mpc. 

7.1 .2. The F2.2 linewidth TF calibration 



The only way to get a match in either case is to vary both the 
gradient and the zero point in the Tully-Fisher relation. We 
find that a good mapping can only be had with TF calibrations 
similar to 



M 



91ogioK-ot - 1- 



But the gradient here, —9, is well outside of the envelope of 
typical i?-band TF gradients and can only be considered as 
extreme. Given that Courteau's 1997 study was primarily de- 



signed to investigate algorithmic definitions of linewidths, and 
that he judged his V2.2 linewidths to be on a par with Hj 
linewidths, it would seem that the only likely explanation for 
the problem lies in the possibility that the sample of chosen 
galaxies is not quiet in the Hubble sense. 



I997| calibration of the Tully-Fisher relation for ^-3- The effect of excluding the edge-on galaxies 



The Courteau 

his V2.2 linewidths, using his H = 70 km/sec/Mpc, is given by 

Mtf = -6.36 (±0.22) logy2.2 - 6.78 . 
We find that the recalibration 

Mtf = -6.55 log ^2. 2 - 7.15 

gives the In A frequency diagram of figure ^ Except for the A 
peak which is non-existent because of lack of data at the dim 
end of the Coruteau sample, we see that the peak structure of 
the hypothesis is, again, reasonably well reproduced. Note that 
the recalibrated gradient is within Courteau's quoted error bar 
This time, the net effect of the recalibration is to increase 
the luminosity of the average galaxy by about 0.77 magni- 
tudes. This is well outside Courteau's quoted error bar for V2.2 
calibrations of 0.46 magnitudes. Since the recaUbrated gradi- 
ent has not changed significantly, this result is consistent with 
the idea that the recalibrated zero point is significantly differ- 
ent from Courteau's original value for the V2.2 calibration. 

7.2. An oddness in the Courteau sample 

As a consistency check on the Vmax calibration, we note that 
the application of this calibration to the Courteau sample with 
H = 70 km/sec/Mpc gives the magnitude mapping 



-23.3,-18.2)^70 



-23.0,-19.7)tf- 



Thus, we see that Courteau's Vmax calibration compresses the 
magnitude range by about 35%, with the dim end being too lu- 
minous by 1.5 whole magnitudes. It is easy to see that varying 
the Tully-Fisher zero point will simply have the effect of shift- 
ing the given TF range (—23.0, — 19.7)tf either up or down - 
but it cannot expand the range to match the Hubble range. 

Similarly, a consistency check on the V2.2 calibration, with 
H = 70 km/sec/Mpc gives the magnitude mapping 



Less that 1% of the Courteau sample consists of edge-on 
galaxies, and so there is no discernible effect arising from their 
exclusion. 
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Fig. 8. The In A distributions for the Courteau 1997 sample 
using his Vmax linewidth definition and H = 70 km/sec/Mpc. 
Vertical dotted lines indicate peak centres of figure 0. Bin 
width = 0.055 



8. The statistical significance of the various 
analyses 

We begin with a broad overview of the various analyses: the 
results of the foregoing analyses, already encapsulated in fig- 
ures g, g, ^ 0, ^ andJU are collected together for convenient 
comparison in table ^The refined specific hypothesis to be 
defined and tested states, briefly, that strong peaks in In A fre- 
quency diagrams of rotation curve data processed in the way 
described, should occur coincidently with the peaks A, B, C 
and D in figure g, which has been derived from our autofolder 
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Courteau sample of 283 objects 



Table 5. In A interval transformation 




Fig. 9. The hi A distributions for the Courteau 1997 sample 
using his V2.2 linewidth definition and H — 70 km/sec/Mpc. 
Vertical dotted lines indicate peak centres of figure 0. Bin 
width = 0.055 

Table 4. Comparison of peak posi- 
tions in In A frequency diagrams for 
MFB, MF, DGHU & SC data 



Sample 


A 


B 


C 


D 


MFB 


3.91 


4.18 


4.73 


5.12 


MFvopt 


3.91 


4.18 


4.68 


5.06 


MFvorig 


3.91 


4.18 


4.68 


— 


DGHU 


— 


4.18 


4.68 


5.06 


oLjYmax 


- 


4.18 


4.73 


5.12 


SCV2.2 


- 


4.18 


4.73* 


5.12 


* indicates 


average 


over three point 


^ 


— indicates 


no significant 


peak 





analysis of Mathewson et al (1992) data. The exact bin-centre 
positions in which these latter peaks lie are given in the first 
row of table ^ whilst remaining five rows give the peak cen- 
tres for the remaining figures. The table makes it clear that the 
peak positions are essentially identical across the four sam- 
ples. In the following, we provide a standardized quantitative 
estimate of the statistical significance of these results. 

8.1. A Standardized methodology 

We noted, in §|l| that an extremely conservative upper-bound 
estimate of the probability of the peaks in figure |l| occurring by 
chan ce alo ne, given the prior hypothes is raise d on the Rubin 
et al 1980 data, was given in Roscoe 1999a| to be less than 
10^^. However, this estimate was derived using a crude ad-hoc 
methodology, and applied specifically to the Persic & Salucci 
folding solution of the Mathewson et al (1992) data. For the 
purpose of enabling statistical comparison across our whole 
analysis, we introduce a standardized methodology - which is 



Interval with 


Interval with 


Associated 


RFT scale 


MFB scale 


peak 


(3.35, 3.65) 


(3.69, 3.93) 


A 


(3.85, 4.15) 


(4.10, 4.34) 


B 


(4.35, 4.65) 


(4.51, 4.75) 


C 


(4.85, 5.15) 


(4.92, 5.16) 


D 



already partly implemented 
each sample. Specifically: 



and apply it where possible to 



- require that all samples to be tested are processed (folding 
etc), as far as is possible, in an identical fashion; 

- define the null hypothesis that all samples are drawn from 
the same background distribution, and that this latter dis- 
tribution is smooth - that is, has no peak structure; 

- set up the specific hypothesis to be tested, and test it via 
a Monte-Carlo simulation which randomly selects a very 
large number of samples from the hypothetical smooth dis- 
tribution. 

The first of these points is, of course, already implemented. 
For the second point, experimentation shows that the final out- 
comes are insensitive to any reasonable choice of 'smooth dis- 
tribution' for the null hypothesis; we define it as the cubic 
spline envelope of the figure distribution shown in figure llOl 



For the third point, the Mathewson et al 1992 sample must be 



tested against the or iginal hypothesis raised on our analysis of 
the Rubin et al 1980 sample, whilst our additional samples will 
be tested against a refined hypothesis raised on the analysis of 
the Mathewson et al 1992 sample. 



8.2. Significance oftlie Matliewson etal \1992 
results: Figure^ 

In essence, we shall revise the original crude ad-hoc estimates 
of the significance of the peaks arising from the Persic & 



Salucci folding of Mathewson et al 1992 data, shown in figure 
in given the original hypothesis raised on the Rubin et al 198C 
data described in fpl. Briefly, this original hypothesis stated 
that, using Rubin scaling, then In A would lie within ±0.15 of 
integer or half-integer values - specifically, the values 3.5, 4.0, 
4.5 and 5.0. The transformation of these Rubin-scaled inter- 



vals to the scaling used by Mathewson et al |1992| is given in 
tableg. The required probability estimate for the peaks of fig- 
ure gwas obtained by generating 10® randomly selected sam- 
ples of 866 measurements e ach (th e number of auto-foldable 
ORCs in the Mathewson et al 1992 sample) from the hypothet- 
ical smooth distribution of figure |10|, and counting how many 
times peaks of the observed (or greater) sizes actually occur in 
the intervals specified in table q[ The frequency at which peaks 
of the required size actually appeared is given in table ^. On 
no occasion did more than two peaks at the required strengths 
appear in the same trial; however, assuming independent prob- 
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Hypothetical smooth distribution 




Fig. 10. Hypothetical smooth In A distribution, defined as a 
cubic spline envelope of actual MFB 1992 distribution. Bin 
width = 0.055 

Table 6. Test of peak-structure in figure g 
Peak frequencies in lO'' trials 



Interval with 


Frequency at 


Peak 


MFB scale 


required 
strength 


label 



(3.69, 


3.93) 


3513 


A 


(4.10, 


4.34) 


758 


B 


(4.51, 


4.75) 


18881 


C 


(4.92, 


5.16) 


66303 


D 



Probability of four sufficiently strong peaks 
appearing in the same trial ~ 3 x 10~^ 



abihties, the observed frequencies allow us to estimate that the 
probability of four peaks of the required strengths appearing 
simultaneously to be estimated at 3 x 10^^. 

8.3. A refined hypothesis for the new samples 

The considerations of §k!3 allow us to refine the hypothesis 
tested there (which arose from a consideration of just 12 Rubin 
et al 



1980 



I objects) as follows: 
The distribution of In A, computed for folded ORCs according 
to the prescription of §pl will show significant peak structure 
with peaks centred on the In A values (3.91, 4.18, 4.73, 5.12), 
where 90% confidence limits on the positions of these peaks, 
calculated in detail in Appendix are given approximately 
by (3.87, 3.98), (4.15, 4.21), (4.70, 4.76) and (5.09, 5.14) re- 
spectively. 



8.4. The significance of the Dale et al results: 
Figure § 



was no freedom to vary this parameter Similarly, the H = 
95 km/sec/Mpc value used to fix the TF zero point arose di 
rectly from optimizing the results of the Mathewson et al 199^ 
analysis given in ffl. Effectively, therefore, the In A frequency 
diagram of figure ffl arose from a single trial and, accordingly, 
the associated probabilities can be calculated directly from the 
figure on the basis of the modified hypothesis of § 5.3. 

The procedure is as before, except that each randomly se- 
lected sample now contains 454 measurements each - which 
is the number of auto-foldable ORCs in the Dale et al sam- 
ple - and the results are shown in table ^ Assuming indepen- 
dent probabilities, the observed frequencies allow us to esti- 
mate that the probability of the three peaks listed appearing 
simultaneously at the observed magnitudes in figure Q to be 
estimated at 8 x 10^^*^. However, it is obvious from table 
that most of the power in this latter result arises from the ex- 
treme naure of peak C. But, even if the number of peak C 
occurrences equalled the number of peak D occurrences, the 
single-trial probablility for this case would still be insignifi- 
cant at 6 x lO"*^. 



8.5. The significance of Mathewson & Forcl \1996 
results: Figure^ 

We now calculate the significance of the peak structure arising 
from our analysis of Mathewson & Ford 1996 data - exhib- 
ited in figure ^ for the_Vppt linewidth results - given the mod- 



The error bars given by Giovanelli et al 1997 for their TF gra- 
dient estimate were so tight (±0.02) that, effectively, there 



ified hypothesis of § |8.3| . The procedure is as before, except 
that each randomly selected sample now contains 1085 mea- 
surements each - which is the nu mber of auto-foldable ORCs 
in the Mathewson & Ford 1996 sample - and the results are 
shown in table 0. 

Assuming independent probabihties, the observed fre- 
quencies allow us to estimate that the single-trial probability 
of four peaks of the required strengths appearing simultane- 
ously to be estimated at 3 x 10^^". However, in this particular 
case: 

- We allowed ourselves the freedom to use an alternative 
linewidth estimation method when the original linewidths 
of Mathewson & Ford 1996 , which gave rise to figure |], 
were inferred to be problematical; 

- Whilst, for the alternative method we used the Giovanelli 
et al gradient value for the TF calibration, we allowed 
ourselves the freedom to search for a new zero point. 
Discounting the fact that we were guided by the magnitude 
mapping technique, it is a fair assessment to say that, typ- 
ically, the search for a suitable zero point required about 
five further independent trials. 

To make a round figure, suppose then that the figure ^ results 
required ten independent trials to obtain. On this basis, we can 
use the single-trial probability and binomial statistics to esti- 
mate the required probability at 3 x 10^^. 
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Table 7. Test of peak-structure probabilities using the refined hypothesis of § 13 





Peaks at observed stren 


gth in lO'' 


simulations 


Single trial 




A 


B 


C 


D 


probability 


Fig 


(Dale) 


- 


11797 


3 


23010 


8 X 10""' 


Fig 


(MF Vopt) 


28106 


131 


1204 


70559 


3 X 10-^° 


Fig 


(MF Vorrg) 


192 


1024 


93262 


- 


2 X 10~** 


Fig 


(SC Vmax) 


- 


44889 


8098 


34653 


1 X lO"'' 


Fig 


(SC 1/2.2) 


- 


44918 


17452 


71252 


6 X lO-'^ 



indicates non-existent peak. 



8.6. The significance ofttie Courteau results: 
Figures^ 

The estimation of the probabilities attached to the Courteau 
results is a little more involved than the other two case. Briefly, 
we argue as follows: 

- We chose at the outset to consider only the Vmax and V2.2 
Courteau linewidth estimates. So there are two degrees of 
freedom arising here; 

- For each of the linewidth estimates, the error bars quoted 
by Courteau for the gradients are approximately ±0.2. 
Given that we allow ourselves the freedom to vary the gra- 
dient within the quoted error bars, this amounts to about 
three independent trials to optimize the gradient value per 
linewidth estimate; 

- Finally, we allowed ourselves the freedom to determine the 
zero points for each of the TF calibrations. Typically, each 
zero point determination required five trials for a given 
gradient determination. 

So, for each linewidth estimate, we can reasonably say that fif- 
teen independent trials were required to determine the optimal 
TF calibration - making thirty independent trials in all. Given 
the single-trial probability of 1 x 10"^ from table It] for the 
Vmax Hnewidth analysis, we can use binomial statistics to es- 
timate the overall probability of obtaining the figure || results 
by chance alone at about 3 x 10"^. 

8. 7. Summary of statistics 

A hypothesis was raised on the basis of a simple analysis of 12 
Rubin et al 198C galaxies; this hypothesis was tested against 
the results obtained from a sample of 866 Mathewson et al 
1992 ORCs and it was confirmed with a probability of « 3 x 



10^ against the obtained results arising purely by chance. 

The analysis of this larger sample of 866 ORCs allowed a 
refinement of the hypothesis, which was subsequently tested 
against the results obtained from three new independent sam- 
plesof 454 ORCs with /-band photometry from Dale et 
al 1997 et seq, 1085 ORCs with /-band photometry from 
Mathewson & Ford |1996| and 283 ORCs with i?-band pho- 
tometry from Courteau 1997 The probabilities of the results 
obtained from each of these samples arising by chance alone 



were estimated at 8 x 10" , 3 x 10"^ and 3 x 10"* respec- 
tively. On the basis of these results, it is reasonable to say that 
the 'discrete dynamical classes' hypothesis for spiral discs has 
been verified at the level of virtual certainty. 



9. A second generation of Tully-Fisher 
methods. 

At a practical level, the foregoing analysis has implications 
concerning the general nature of Tully-Fisher methods which 
we shall expand upon in this section. These can be summa- 
rized as: 

1 . The discrete dynamical states phenomonology provides an 
absolute standard whereby Tully-Fisher methods can be 
absolutely calibrated on a large enough sample; 

2. The Tully-Fisher method is shown to have two equivalent 
formulations, one of which is the familar one; 

3. The Tully-Fisher relation is shown to be augmented by a 
corresponding relation which defines where on an ORC a 
linewidth should be measured. Thus, a means is provided 
whereby Tully-Fisher calibrations and linewidth determi- 
nations can, in principle, be algorithmically, absolutely 
and simultaneously determined from any given large sam- 
ple of ORCs; 

4. The concept of 'linewidth' is shown to be not absolute - 
that is, linewidth can be defined with virtually unlimited 
freedom. But, once defined, the form of the coiTesponding 
Tully-Fisher calibration becomes fixed. 

9. 1. Thie essential background 

The consequences listed above rest upon the assumption that 
the power-law analysis can be shown to provide a high-quality 
resolution of ORC data in that part of the disc where the dis- 
crete dynamical classes phenomonology manifests itself. To 
demonstrate this, the first thing of significance to be consid- 
ered is the {a, In A) pl ot, given in figure [ill for the combined 
Mathewson et al |1992| , Mathewson & Ford |1996| , Dale et al 



1997 et seq samples of 2405 objects having /-band photom- 
etry and reliably foldable ORCs. We now make a detailed 
analysis of the 1951 (a, In A) points in figure pi] which are as- 
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MFB + MF + DGHU samples 



5 - 
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Fig. 11. Plot of (In A, a) for 2405 galaxies 



Tables. In A model 



Variable 


Coeff 


t ratio 


const 


-1.596 


-11 


M 


-.316 


-46 


S 


-0.0002 


-2 


a 


7.614 


28 


aM 


0.474 


34 


aS 


0.0050 


18 



N=1951;7?^ = 92% 



sociated with the specific subset of the Mathewson et al 1992 
+ Mathewson & Ford 1995 ORCs. In the context of the model 



V_ 

we find 

In A = 



Rf) 



A = 



ID a 



(10) 



InVo — a In Rq 
-1.596 -0.316 A/ 
+0.0050 a S" 



7.614 a + 0.474 a M 



(11) 



for the two /-band Mathewson samples. Here M is absolute 
magnitude, and S is surface brightness defined as average so- 
lar luminosities per square parsec for the whole area inside 
the optical radius, R^s, as defined by, for example, Persic & 
Salucci 1995 This particular model was obtained by using the 



original Mathewson linewidths for the Mathewson et al 1992 
sample, and the Vopt linewidths for the Mathewson & Ford 



1996 sample, and by rejecting all 3 a outHers (about 4% of 



the total), and it accounts for about 92% of the variation in 
figure O. The detailed statistics are given in table K^ and it is 
clear that, except for surface brightness, S, all of the chosen 
predictors have an extremely powerful effect in the model. 



9.2. Extraction of models for In Rq and InVo 

Introducing the 'discrete dynamical states' phenomonology 
into the discussion, we now note that the model can be decom- 
posed into separate expressions for (In Ro,\nVo) according to 



\nAi = lnVo-alni?o, i = 1,2,3,4 
\nVa = -1.596 - 0.316 Af + r 



(12) 



lni?o = -7.614 - 0.474 A'f 



0.0050 5*+- 
a 



where In Ai = 3.91, 4.18, 4.73, 5.12 and denotes the posi- 
tions of the A,B,C,D peaks, and F is an arbitrary function 
(or parameter) arising from the decomposition of In A;. It is 
easily verified that ( |lO| ) is invariant with respect to an arbitrary 
choice of this F. 

An effective visual verification of the fit of model dl^) is 
obtained by noting that if we use the definitions of ( p^ in 
the dimensionless form of (pX3|), and then regress lii(T^/Vb) 
on \n{R/Ro), then we should find a null zero point for each 
ORC - except for statistical scatter. Because of the invari- 
ance of ( pX)[ ) with respect to F, the results of this regression 
will be independent of any chosen F, and so we can set it to 
zero for this particular exercise. Figure [Tj (left) gives the fre- 
quency diagram for the actual zero points computed for the 
combined Mathewson et al (1992, 1995) samples from which 
the model (O) was derived, whilst a wholly independent test 
of the model is given by figure O (right) which gives the fre- 
quency diagram for the zero points derived from the Dale et 
al 1997 et seq sample using the model (pi]). It is clear that 
there is absolutely no evidence to support the idea that these 
zero points are significantly different from the null position - 
so that ( p^ with ( 12 ) can be considered to give a very effective 
resolution of ORC data in the disc. 



MFB+MF zero points 
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Fig. 12. Plot of zero point for 1951 MFBh-MF galaxies and 
454 DGHU galaxies 
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9.3. Implications 1: The classical Tully-Flsher 
relation 

9.3.1 . Case r = 0: Retrieval of the standard 
formulation 

We begin by noting that a rearrangement of the second equa- 
tion of (|l2|), converting In Vq to log Vq, and setting F = 
gives 



M = -7.29 log Vb - 5.05. 



(13) 



Interpreting the characteristic scaling velocity Vq as the disc 
'rotation velocity', we can immediately recognize this rela- 
tion as similar to a typical /-band Tully-Fisher calibration - 
but with a gradient on the low side. However, this latter re- 
lation was derived by combining two samples having differ- 
ent linewidth characteristics to make a specific point about the 
general analysis. Consequently, we might expect to find more 
typical gradients by analysing the four samples individually. 
In fact, we find: 

Mi = -7.72 log t/o - 3.79 

Mi = -7.85 logK) - 3.54 

Mi ^ -7.86 log K) - 3.71 

Mr = -6.56 log K) - 6.48 



respectively for the Mathewson et al 1992 , Mathewson & Ford 
1996|, Dale et al |1997| and Courteau]l997| samples - the latter 



using Courteau's Vmax linewidth definition as an example. We 
see that the three /-band gradients are now perfectly typical /- 
band TF gradients, whilst the single /?-band gradient is now a 
perfectly typical /?-band TF gradient. Similarly, for each case, 
the zero points are typical for H w 85km/sec/Mpc. We can 
therefore reasonably conclude that the second equation of ( [l2| ) 
with F = recovers the standard formulation of the Tully- 
Fisher relation. 

It now follows directly that the third equation of (|l2|), for 
the characteristic scaling radius Rq, effectively defines the po- 
sition on an ORC at which the rotation velocity, Vq, is to be 
measured. That is, for the example given, the F = case 
corresponds to measuring a characteristic rotation velocity 
In Vo = —1.596 — 0.316 M at a characteristic radial position 
In Ro = -7.614 - 0.474 M - 0.0050 S. 

9.3.2. Case F / 0: Arbitrariness in calibration 
procedures for TF 

Given the F = case, it is now easy to see that the F ^ 
case corresponds to measuring a characteristic rotation veloc- 
ity, V = exp(F) Vq, at the radial position R = cxp(r/a) Rq. 
Thus, according to the present considerations, there is an ar- 
bitrariness in calibrating Tully-Fisher relations which consists 
in the freedom to choose which point on an ORC to use as the 
characteristic radius for the measurement of the characteristic 



rotation velocity. Also, since F is an arbitrary function - and 
not just a constant - this freedom of choice can extend to, for 
example, a procedure whereby the place at which the rotation 
velocity is defined can be luminosity dependent. 

This has significant consequences since, as mentioned in 



B.2, there are several non-equivalent heuristically defined 



methods which authors use to estimate optical linewidths and 
there is, as yet, no concensus about which method is 'best'. 
But, according to the above comments, it is clear that the 'best' 
such method is the one which most consistently selects the 
radial positions R = exp(r/a) Rq across the sample under 
consideration, for some F. 

9.4. Implications 2: Non-classical Tully-Fisher 
formulations 

The foregoing considerations give rise to distinct Tully-Fisher 
formulations neither of which is logically prior to the other. 
They are equivalent in the sense of being equally applicable to 
the problem of scaling. 

9.4.1. The quasi-classical formulation 

If we choose the F = calibration, then (n2b can be expressed 
as: 



InK) = -1.596 -0.316 A/ 

ln/?o = -(1.596 + 0.316 M + In A,)/a 
i ^ 1,2,3,4. 



(14) 



As we saw in § 9 . 3 . 1 , the first of these equations provides a typ- 



ical /-band Tully-Fisher calibration whilst the second effec- 
tively says where on an ORC the linewidth, Vq, is to be mea- 
sured. Note that, in this formulation. In Vq varies only with M 
as in any conventional Tully-Fisher calibration, but Rq varies 
with M and the dynamically measured quantities (a, In A). 

9.4.2. The non-classical formulation 

Equivalently, with the same F = calibration, ( [l2[ ) can be 
expressed as 

InVb ^ -a(7.338 +0.461 M + 0.00449 5) + In A, 

i = 1,2,3,4 
ln/?o = -(7.338 +0.461 M + 0.00449 5). 

In this formulation, we see that the classical situation is re- 
versed. Now In Vq varies with M and the dynamically mea- 
sured quantities {a, In A). Consequently, unlike in the classi- 
cal formulation, both gradient and zero point will vary across 
the luminosity/dynamical range of a sample. By contrast, it is 
Rq which varies with just luminosity properties, M and S. 
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9.5. Comments on general calibration procedures 

For the sake of explicitness, we shall use ( [l4| ) to outline a sug- 
gested algorithmic approach to the calibration procedure for a 
sufficiently large sample. The approach is non-trivial and has 
yet to be developed into a fully working method: 

1 . Use the model (|l4|), or something similar, as first guess at 
the final calibration; 

2. Fold an unsealed ORC and determine its a - which is a 
scale-independent quantity; 

3. Estimate which In Ai is appropriate for the ORC using the 
(a : In A) correlations implicit in figure [ill; 

4. Iterate on the two equations of ( p4} ) and the power-law 
V = AR" fitted to the ORC, to find {Rq, Vq, M) for the 
galaxy concerned. The individual ORC can now be scaled; 

5. Complete steps 2,3 and 4 for the whole sample and then: 

(a) calculate the new In A frequency diagram and compare 
with some template distribution such as that of figure 

i 

(b) check if the calibration which comes out of the model 

(in the manner of (Ulb) is the same as the model which 

went in at step 1 - to within tolerance; 
If both these tests are satisfied, then accept the current cal- 
ibration. If not, adjust the calibration, and repeat steps 2, 
3, 4 and 5 until converged. 

10. Major Stability issues 

Apart from the importance of good Tully-Fisher calibrations 
and linewidth determinations, the successful extraction of the 
discrete dynamical states phenomonology from ORC data is 
also critically dependent on the quality of the folding process 
and the algorithmic computation of In A. We consider these in 
turn. 



MFB sample with MFB folding solution 




Fig. 13. In A distribution for the Mathewson et al 1992 sample 
with MFB folding and original Mathewson et al TF scaling; 
Vertical dotted lines indicate peak centres of Persic & Salucci 
solution. Bin width = 0.055 

the effect of discarding, typically, the least accurate 40% of 
velocity determinations on an given ORC. They then applied 
a labour-intensive 'by-eye' approach to the folding problem. 
The folding algorithm developed for the present study (Roscoe 



1999c ) followed Persic & Salucci in utiUzing both of these re- 



finements. 

To illustrate the effect of folding without these refine- 
ments, figure im shows the In A frequency diagram computed 
for the Mathewson et al sample 1992, but using Mathewson 
et al's own folding solution. This is to be compared directly 
with figures |l] and Q. It is quite clear that the discrete states 
phenomonology is heavily obscured in figure [l^. 



10.1. The quality of folding 

In the literature, the folding of ORCs is most commonly de- 
scribed in the context of estimating linewidths for TF appli- 
cations and, in such circumstances, the folding process gen- 
erally amounts to determining an estimate of Vsys, the ve- 
locity equivalent of the syste matic redshift. So far as we are 
aware, Persic & Salucci 1995 were the first to address the idea 



that ORCs folded for TF purposes were not necessarily folded 
with sufficient accuracy for the purpose of studying the inte- 
rior dynamics of galaxy discs. It was for this reason that they 
provided their own folding solutions for the Mathewson et al 



1992 sample. Possibly their main innovations to the folding 
process were, firstly relaxing the assumption that the dynamic 
centre of a spiral coincided with its optical centre. Hence, their 
folding process became a two-parameter problem - to deter- 
mine Vsys and the angular offset between the dynamical and 
optical centres, A$ say. Secondly, they introduced the pre- 
folding data filtering technique described in §3.1 which has 



10.2. The quality of In A determinations 

The quality of In A determinations is equally critical to the 
extraction of the discrete states phenomonology. Specifically, 
I n A is computed according to the 'black-box' algorithm of 
§3.3 which has the effect of, on average, discarding approxi- 
mately the inner 10% of any given ORC. The value of In A is 
then computed directly on the remaining outer segment of the 
ORC. 

The typical consequences of not using this algorithm, and 
computing In A on the whole ORC, are shown in figure nj 
which gives the In A frequency diagram arising from the 



Mathewson et al 1992 sample after it has been folded using 



our autofolder technique but without using the computational 



method of 5:3.3. Again, we see that the discrete states phe- 



nomonology is completely obscured. 
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MFB sample without Ln(A) algorithm 
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Fig. 14. In A distribution for the Mathewson et al 1992 sample 



with auto-folding and original Mathewson et al TF scaling, but 
without using the In A algorithm; Vertical dotted lines indicate 
peak centres of Persic & Salucci solution. Bin width = 0.055 

10.3. Effect of folding refinements on zero-point 
determinations for scaled ORCs 

It is of interest to consider the zero-point distribution arising 
from regressing the scaled velocities, ln(V^/Vb) on the scaled 
radial coordinates, h\{R/Ro), when the {Rq, Vq) model is ob- 
tained from ORCs folded using simply Vsys determinations 



and without applying the In A algorithm of §3.3. To illus- 
trate the point, the right panel of figure |l5| show the distri- 
bution arising from the combined Mathewson et al 1992 and 
Mathewson & Ford 1996 samples using their original folding 
solutions, and without using the hi A algorithm. The left panel, 
included for comparison, arises from the same combined sam- 
ple, but with our various folding refinements included. 

It is obvious from the figure that the effect of using the 
original folding solutions and omitting our refinements is sim- 
ply to add unbiased noise to the zero-point solutions. However, 
it is interesting to note that, even when the original folding so- 
lutions are used and our various refinements are omitted, the 
zero-point frequency distribution still provides strong support 
for the power-law description of ORCs. 



11. General dynamical Implications 

The present analysis has allowed us to established to a very 
high degree of statistical certainty that the parameter h\A 
appears constrained to take on discrete values ki,k2,---- 
Combining this with, say, ( pi] ) which showed In A w 
F{M, S, a), then we must have F{M, S, a) = ki,k2, .... 
Thus it appears that spiral galaxies are constrained to exist on 
one of a set of discrete class planes in the three-dimensional 
(M, S, a) space. This then gives rise to one of two broad pos- 



Fig. 15. Left panel: zero-point distribution arising from 
mfb-nmf combined sample of scaled ORCs with auto-folder 
and In A algorithm. Right panel: similar to left panel, but with 
Mathewson's original folding solution and no In A algorithm. 
Bin width = 0.055 



sibilities: at some stage in its evolution a spiral galaxy some- 
how moves onto one of these discrete class planes and then: 

- is either necessarily constrained to remain on this plane 
over the whole of its evolution; 

- or has the possibility of transiting to other planes in very 
short periods of time, so that the planes themselves repre- 
sent an evolutionary sequence. 

In either case, the primary difficulty is identifying mechanisms 
which might generate such discrete sets of possible dynamical 
class planes and, in the following, we make three conjectures. 

11.1. First conjecture 

In the wider world on non-linear dynamical systems, it is not 
difficult to find (or to construct) systems for which consistent 
solutions can only be found when certain algebraic conditions 
are satisfied - not eigenvalue problems, but analogous to them. 
If, in such a system, the algebraic condition had the form of 
a quartic defined in the system's parameter space then, im- 
mediately, one would have a situation in which the solutions 
fell into four discrete dynamical classes. Our first conjecture 
is, therefore, that disc dynamics are governed by just such a 
non-linear system so that the discrete dynamical class phe- 
nomonology is merely a manifestation of an unknown alge- 
braic consistency condition. 



11.2. Second conjecture 

Within the context of inflationary cosmologies, scalar fields 
play a critical role in the early universe. Since it is not dif- 
ficult to induce oscillatory behaviour in these scalar fields, it 



18 



D. F. Roscoe: Discrete Dynamical Classes For Galaxy Discs 



becomes natural to consider the possibility that the four dy- 
namical classes of our analysis are the frozen imprint of four 
distinct galaxy-formation epochs in the early universe, each 
associated with a distinct oscillation of the scalar field. If this 
is the case, then the phenomonology provides a very strong 
constraint on these early galaxy formation processes. 

11.3. Third conjecture 

Until recently, the most popular candidate for constituting the 
dark matter that is required to produce the generality of ob- 
served rotation curve shapes (especially for the dimmer galax- 
ies, and beyond optical limits in general) has been the neu- 
tralino. However, this conjecture is now in serious difficulties 
since simulations show that these massive particles would tend 
to accumulate in large agglomerations which would cause se- 
rious disruption of galactic structure which is not observed. In 



response to this state of affairs, Arbey et al 2001 have proposed 



the existence of a massive non-self interacting charged scalar 
field as an alternative source of 'dark-matter' action. They pro- 
vide a fairly comprehensive analysis to show how their conjec- 
ture provides extremely good fits to the six brightest of Persic 



et al's 1996 eleven classes of universal rotation curves. 



More particularly, from our point of view, the conjectured 
dark halo of Arbey et al 2001 is based on a primitive 'boson 
star' model, and this model is constrained so that it gives rise 
to a sequence of discrete energy eigenstates, each of which 
is associated with distinct measures of rotation for the 'boson 
star' . Whilst it is true that these discrete rotation states do not 
correspond in any direct way to the discrete dynamical classes 
discussed in this paper, it must be considered remarkable that 
such a prediction arises at all in a paper which was primarily 
directed towards attempting to model the universal ORCs of 



Persic et al 1996. It seems, therefore, that the phenomonology 



discussed here may potentially be understood in terms of the 
quitessential galactic haloe. 



12. Summary and Conclusions 

We have analysed four separate large ORC samples to show 
that the discrete dynamical classes hypothesis for galaxy discs 
is supported by the data at the level of virtual certainty. The im- 
mediate significance of the phenomonology is that any given 
spiral galaxy appears to be constrained to evolve over one of a 
discretely defined set of dynamical class planes, existing in a 
three-dimensional (M, S, a) space where Al is absolute mag- 
nitude, S is surface brightness and a is a parameter computed 
for each galaxy from its rotation curve. 

We have shown how the broader analysis implies the exis- 
tence of a second generation of Tully-Fisher methods which 
augment the classical Tully-Fisher relationship with a sec- 
ond, similar, relationship which defines where on an ORC the 
linewidth is to be measured. An algorithm for the absolute cal- 



ibration of this extended Tully-Fisher method was then pro- 
posed. 

We have conjectured three possible mechanism for 
this phenomonology, one based on the idea that the phe- 
nomonology is the manifestation of an algebraic consistency 
condition in some (unkown) non-linear dynamical system the- 
ory which describes disc dynamics, a second based on the no- 
tion of a sequence of distinct galaxy-formation epochs in an 
oscillating early universe, and a third based on the dynami- 
cal effects of quintessential halos acting as the source of dark 
matter around spirals. 

Whatever the truth of the matter, it seems certain that the 
existence of the distinct dynamical classes poses very difficult 
questions for the standard galaxy formation theories, and will 
have a a potentially profound affect on our developing under- 
standing of galactic dynamics and evolution in particular, and 
the cosmos in general. 
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Appendix A: Tlie Effects of Tully-Fisher Scatter 
on Inyl Profiles 

All the distance scaling in this analysis is performed using 
Tully-Fisher methods, which possess well-understood inher- 
ent sources of error Consequently, we need to understand the 
extent to which these errors can affect the phenomenon which 
is the subject of the present analysis. 

Mathewson et al 1992 report a magnitude scatter of about 
0.32 for their sample (our best) which compares favourably 
with that of 0.35 reported by Courteau for his sample. These 
correspond to a scatter of less than 20% for distance measure- 
ments and, in the following, we analyse the effects of such un- 
certainties on our proposed analysis to demonstrate that they 
cannot wash out any potential peak structures of the type seen 
in figure nl 

Suppose that each galaxy in the sample has had its distance 
exactly determined, and that R denotes the corresponding ex- 
act radial scale. Then V = AR" implies 

InV = \nA + a\nR. 

The existence of uncertainties in the Tully-Fisher distance 
scale with a typical scatter of 20% can be accounted for by 
the replacement R -^ kR where 0.8 < k < 1.2, so that 

InV = \nA + a\nk + alnR = lnA' + a\nR, 
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where In^' = \nA + a\nk. We immediately see that un- 
certainties in the distance scale affect the zero point in the 
(In R, In V) relationship, but leave the gradient a unaffected. 
Since 0.8 < fc < 1.2 then, as an approximation, —0.2 < 
hifc < 0.2 so that \nA' = liiA± 0.2a. The peak structures 
of figure n^ lie in the approximate range 3.9 < In A < 5.1 



and (b) algortihmically determined linewidths based on 
Vopt, the velocity at the optical radius as defined by Persic 



(that is, 1.7 < log^ < 2.2) and reference to Roscoe 1999b 



(figure 8) shows that this corresponds to the approximate range 
0.18 < a < 0.55, where low a corresponds to the brightest 
galaxies and vice versa. Therefore, for the brightest galaxies, 
we have \nA' = In A ± 0.04 whilst, for the dimmest galaxies, 
we have \nA' = \nA± 0.1. That is, uncertainties of 20% in 
the Tully-Fisher distance scale create uncertainties in In A of 
±0.04 at the bright end, and of ±0.1 at the dim end. 

It follows that, since the mean separation of the In A peaks 
in figure [l| is about 0.4, these uncertainties in the Tully-Fisher 
distance scale are incapable of washing out the discrete peak 
structure observed in the In A distribution. This analysis pro- 
vides the confidence required to analysise further samples on 
the same basis. 

Appendix B: The elimination of artifact as a 
meclianism 

It is necessary to be as certain as is possible that the phe- 
nomonology being claimed is not created as an artifact of any 
particular procedure or data sample. There are three indepen- 
dent routes by which such an artifact (no matter how remote 
the possibility) could infiltrate the process, and these can be 
Usted as: 

- the original process of measuring ORCs; 

- the method of linewidth estimation; 

- the folding process. 

6. 1. The original measuring process 

The possibility that an artifact can enter via the ORC measure- 
ment process is minimized by the fact that we analyse four 
different samples originating with three distinct groups of as- 
tronomers using three different telescopes in different hemi- 
spheres. A detailed discussion of the samples is given in §0. 

B.2. Linewidtli determinations 

The possibility that an artifact can enter via the linewidth es- 
timation process is minimized by the fact that five different 
methods of linewidth determination are employed over our 
four samples: 



- The Mathewson et al 1992 sample has linewidths de- 
termined using a non-algorithmic and intuitive 'eye-ball' 
method (private communication); 

- the Mathewson & Ford 199f sample is considered using 
(a) Mathewson's own 'eye-balF linewidth determinations 



ropt 

& Salucci 



1995; 



- the Dale et al 



1997 



sample has linewidths determined by 
an elaborate process based on the Vqo% — V^io% linewidth 
estimator, which is similar to the technique first introduced 
by Dressier & Faber|1990|; 



- the Courteau sample was taken from the study of Courteau 



1997 which was explicitly designed as a study of objec- 



tive black-box methods of optical linewidth estimation. He 
tests a variety of methods, and we present results using 
those two which he judges to be the best, V2.2, and the 
worst, Vmax, respectively. 

B.3. Tlie folding process 

Finally, the possibility that an artifact can enter via the fold- 
ing process is minimized by the fact that the phenomenon is 
observed when either of two quite distinct folding methods is 
used. These two methods, described below, share the features 
introduced by Persic & Salucci 1995 that (a) the folding pro- 
cess requires the determination of two parameters, the primary 
one being Vgys the systematic velocity (which is usually the 
only one considered), and the secondary one being A$ which 
measures the angular offset between the optical and the kine- 
mati c centres and (b) the pre-folding data filter described in 
§3.1 These features are considered necessary for the purpose 
of maximising folding accuracy. 

B.3.1 . The method of Persic and Salucci 



Persic & Salucci |1995| were primarily interested in using ro- 
tation curves for studies of the interior dynamics of spiral 
galaxies and so, by their own criteria, had a requirement for 
a large sample of particularly accurate ly fol ded ORCs. They 
took the 965 ORCs of Mathewson et al 1992| and used an eye- 
ball method of folding to produce a sample of 900 good-to- 
excellent quality folded ORCs; as a qualitative measure of the 
effort expended to produce this sample, we can note that it 
took these two authors about a year to process it (private com- 
munication). Every velocity measurement in the Mathewson 
et al 1992 sample came provided with a parameter (varying on 
the range (0, 1)) which estimated the relative internal accu- 
racy associated with the measurement. Persic & Salucci 1995 
found that the accurate folding of any given ORC required the 
rejection of any individual velocity measurement for which the 
associated accuracy parameter was < 0.35. In the present con- 
text, only the Mathewson et al 1992| sample has been folded 
with this method. 



B.3.2. The auto-folder method of Roscoe 1999c 



This method was developed in anticipation of the need accu- 
rately to fold the Mathewson & Ford 1996 sample of 1200h- 
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ORCs on a reasonable time-scale. The details of this method 

but, briefly, it is based on the 



1999c 



are described in Roscoe 
formal minimization of the symmetric components in Fourier 
representations of ORCs with respect to variations in the two 
folding parameters. 

The folding method was developed on that subset of the 
Mathewson et al 1992 sample used by Persic & Salucci 1995 



and, corresponding to the experience of these latter authors, 
we found that the optimal trade-off point between the qual- 
ity of individual velocity measurements, and the volume of 
good-quality data available for the automatic folding method, 
required the prior rejection of any individual velocity mea- 
surement which had an associated relative accuracy parameter 

< 0.4. This is roughly equivalent to the requirement that the 
absolute error of any given velocity measurement should be 

< 5%. 

This folding metho d has been used here on the samples of 
Mathewson et al 



1992, Mathewson & Ford 1996 



19971 et seq and Courteau [1997 



Dale et al 



Appendix C: Overview of thie ORG dynamicai 
partitioning process 

We discuss two methods of assessing the efficiency and effec- 
tiveness of the process described in §3.3 for the computation 
of In A and establish, with virtual certainty, the truth of the 
statements that: 

- the innermost parts of ORCs exhibit behaviour which is 
qualitatively sharply distinguished from that exhibited by 
the outermost parts of ORCs; 

- the size of such innermost sections can be quantified in 
terms of a radial measure, Rmin say, which can be shown 
to be extremely powerfully correlated with the indepen- 
dently defined optica l radius, Ropt (here, as given by 
Persic & Salucci 1995 of the disc. Since Ropt carries phys- 



ical information about the system, then we must conclude 
that the algorithmically estimated Rmin Hkewise carries 
physical information about the system. Given the qual- 
ity of the statistics involved, these two points are entirely 
sufficient to establish that Rmin does, in fact, define a 
real boundary between distinct dynamical regimes which, 
in turn, gives a concrete justification to the technique by 
which it is estimated. We suggest that Rmin could act as a 
tracer for the gravitational radius of the core on the basis of 
the circumstance that there appears to be no other possible 
interpretation. 



C. 1. The first test of dynamical partitioning 



We argued in Roscoe 1999b| that, if Rmin really was a tracer 
for the gravitational radius of the core, then we might expect 
to find a positive correlation between Rmin and Ropt - on the 
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Fig. C.l. Scatter plot of Ropt : Rmin for MFB, MF DGHU 
and SC data with dynamical partitioning 

grounds that galaxies with large cores might be expected to 

have large optical radii etc. 

The top-left panel of figure ZA_ sh ows th e (Ropt ■ Rmin) 
scatter plot for Mathewson et al data 1992 and we note the 
existence of an extremely strong positive (Ropt ■ Rmin) cor- 
relation. The top-right pane l shows an equally strong effect 
for Mathewson & Ford 1996 data, whilst the botto m-left panel 
shows a similarly strong effect for the Dale et al 1997 et seq 



sample, and and the bottom-right panel shows the same effect 
for the much smaller Courteau sample. The statistics underly- 
ing these diagrams are given in the various columns of table 
C.l and these confirm in quantitative terms what is obvious in 



the diagrams. We can therefore conclude with certainty that, 
since Ropt is an objectively defined physical boundary, the 
Rmin is similarly an objectively defined physical boundary. 
We speculate that, in fact, it is a noisy tracer of the boundary 
between core-dominated and disc-dominated dynamics in the 
disc. The foregoing considerations lead to the following con- 
clusions: 

- The application of the dynamical partitioning process pro- 
duces a very powerful, although noisy, Rmin '■ Ropt cor- 
relation confirming that Rmin (as computed by dynamical 
partitioning) is a strong tracer for Ropt', 

- The computed value of Rmin defines a physical tran- 
sition boundary between core-dominated dynamics and 
disc-dominated dynamics. 

Taking these items together, and noting the absence of any 
other obvious interpretation, we conclude that Rmin almost 
certainly represents a dynamically derived tracer of the gravi- 
tational radius of the core. 



C.2. Tfie second test of dynamical partitioning 

The second definitive formal measure of the statistical effi- 
ciency of dynamical partitioning is given in table |C.2|. The first 
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Table C.l. 


Ropt — bo + bi 


tirnin 








MFB MF 


DGHU 


SC 


Predictor 


Coeff t-ratio 


Coeff t-ratio 


Coeff t-ratio 


Coeff t-ratio 


Const. 


5.26 44 


5.95 55 


5.97 23 


12.45 25 


^min 


1.74 28 


1.54 30 


1.06 11 


1.21 12 




R^ = 48.9% R^ = 45.7% 


R^ = 22.2% 


i?" = 35.5% 



Table C.2. Effect of dynamical partitioning technique on power-law fits 



Dynamical 
partitioning 


MFB 

Mean rms 


MF 
Mean rms 


DGHU 

Mean rms 




SC 

Mean rms 


Before 
After 


10.8 X 10"^ 
2.4 X 10"^ 


12.0 X 10"^ 
3.4 X 10"^ 


1.1 X 10"^ 
0.6 X 10"^ 




1.8 X 10"^ 
1.4 X 10"^ 


Statistics calculated on the 866 ORCs of MFB, the 1085 ORCs 
the 454 ORCs of DGHU and the 283 ORCs of SC. 


ofMF 





row of this table gives the averaged root mean square (rms) 
eiTor calculated from fitting power-laws to each of the 866 



foldable ORCs of the Mathewson et al |1992 samp le, the 1085 

1996| sample, the 



foldable ORCs of the Mathewson & Ford 



454 foldable ORCs of the Dale et al |1997|et s eq sample and 
the 283 foldable ORCs of the Courteau |1997| sample before 
the dynamical partitioning process. The second row gives the 
corresponding averaged rms values after the dynamical par- 
titioning process. Column 2 refers to Mathewson et al 1992 
data, column 3 refers to Mathewson & Ford 1996| data, coul 
umn 4 refers to Dale et al 1997 et seq data and column 5 refers 
to Courteau 1997 data. The Mathewson et al 1992 data shows 



an almost 80% reduction in its mean rms, the Mathewson & 
Ford 1996| data shows a 72% reduction in its mean rms, the 
Dale et al 1997 et seq data shows a 45% reduction in its mean 
rms and the Courteau 1997 data shows a 29% reduction in its 
mean rms. 

Since the dynamical partitioning process leads to the dis- 
carding of only 12% of the Mathewson et al 1992 and the 
Mathewson & Ford |1996| data (out of a total of about 37000 



individual measurements over the two samples), of only 9% 
of the Dale et al 119971 et seq data (out of a total of 15000 indi- 



vidual measurements) and of only 10% of the Courteau 1997 



data (out of a total of about 17000 individual measurements), 
then we can categorically state that the table provides con- 
clusive evidence for the statement that the deviation from the 
power-law fit is strongly concentrated on the inner 10% or so 
of ORCs. From this we can conclude that the behaviour of the 
inner 10% or so of ORCs is qualitatively sharply distinguished 
from the outer 90% or so, as we would expect on purely dy- 
namical grounds. 



Appendix D: Confidence iimits on peai< 
positions 

We use bootstrap techniques to show that the uncertainties in 
the positions of the A, B, C, D peaks of figure 0, are small. 
The process adopted is described as follows: 

- Partition the range of In A into the sub-ranges 
(2.200,4.070) containing peak A, (4.070,4.455) 
containing peak B, (4.455, 4.950) containing peak C and 
(4.950, 6.000) containing peak D; 

- Use the actual In A data-set of N distinct In A values to 
generate 1000 bootstrapped simulated data-sets, each con- 
sisting of N values generated by random selection from 
the real data-set with replacement; 

- Form the frequency diagram for each of the 1000 simu- 
lated data-sets, and record the position of the largest signal 
only in each of the four sub-ranges defined above; 

- Form the frequency diagram for this 'largest signal' data 
set - this is plotted in figure p.l[ 

We now discuss the peaks in figure p.l| in turn: 
PeakA:lnA = 3.91 

This peak coincides exactly with peak A of figure 0. Since In A 
correlates in a strong positive sense with absolute luminosity 



(Roscoe 1999b then this peak corresponds to the dimmest end 



of the sample, so that greater measurement uncertainties prob- 
ably account for the relative broadness of this peak; even so, 
this peak is still very tightly defined and the 90% confidence 
interval this peak is (3.87,3.98). 
Peak B:lnA = 4.18 

This peak coincides exactly with peak B of figure g and is 
the most straightforward case, being very tightly defined with 
no subsidiary peaklets; in fact, more than 97% of the sampled 
B peaks fall in a single bin, so that we can assert that an ex- 
tremely conservative 90% confidence interval for this peak is 
given by the bin boundaries as (4.15,4.21). 
Peak C:lnyl = 4.73 
This peak exhibits a bi-modal structure with the major mode 
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MFB Bootstrap frequency diagram 
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Fig.D.l. Frequency diagram for tlie positions of the largest 
signals near the original A, B, C, D peaks, bootstrapped from 
the MFB sample over 1000 trials 

coinciding exactly with the peak C of figure 0; the minor peak 
contains about 10% of the sampled C peaks, with the remain- 
der being in the major peak. The boundaries of the major peak, 
given by (4.70,4.76), therefore give an estimate of a 90% con- 
fidence interval for the C peak. 
Peak£»:lnyl = 5.12 

Again there is a bi-model structure, but this time the minor 
mode is very much more prominent than it is for the C peak. 
We find that the minor peak contains approximately 20% of 
the sampled D peaks, with the remainder being in the major 
peak. To simplify the complexities presented by this bi-modal 
structure, we simply choose the boundaries of the major peak 
as estimates of the required 90% confidence interval for the D 
peak. 
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